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Abstract 

Using the non-relativistic hydrodynamic limit, we solve equations of motion 
for Einstein gravity and Gauss-Bonnet gravity with a negative cosmological 
constant within the region between a finite cutoff surface and a black brane 
horizon, up to second order of the non-relativistic hydrodynamic expansion 
parameter. Through the Brown- York tensor, we calculate the stress energy 
tensor of dual fluids living on the cutoff surface. With the black brane solutions, 
we show that for both Einstein gravity and Gauss-Bonnet gravity, the ratio of 
shear viscosity to entropy density of dual fluid does not run with the cutoff 
surface. The incompressible Navier-Stokes equations are also obtained in both 
cases. 



1 Introduction 

The AdS/CFT correspondence [1, 2, 3, 4] relates gravity in an anti-de Sitter (AdS) space- 
time to a strongly coupled conformal field theory (CFT) living on the boundary of the AdS 
space. Recently the AdS/CFT correspondence has been applied to various fields. By use 
of the AdS/CFT correspondence, one can calculate many quantities of strongly coupled 
CFTs through dual gravity theories. A remarkable example is the calculation of the ratio 
of shear viscosity to entropy density r]/s of some field theories dual to the AdS Einstein 
gravity [5, 6, 7, 8]. In Einstein gravity, it was found that the ratio is a universal value 
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1/4:71, while in the case with R 2 corrections [9, 10, 11, 12] there is a negative additional 
correction term. Furthermore, it was shown that rj/s only depends on the value of effec- 
tive coupling of transverse gravitons evaluated on the horizon [13, 14, 15, 16]. Under in 
the hydrodynamic limit, it was generally proven that some linear response coefficients are 
universal both in the AdS/CFT correspondence and membrane paradigm [14]. In partic- 
ular, one can consider a fictitious membrane at a constant radial radius to express the 
AdS/CFT response in terms of the membrane paradigm language [14]. The dependence 
of the diffusion constant of dual fluid on the cutoff surface is interpreted as the Wilson 
renormalization group flow [17]. Some studies relating the radial radius of AdS space to 
energy scale of dual CFT appear in [18, 19, 20]. Some forms of holographic renormalizaton 
group flow equation independent of the cutoff were given recently in [21, 22, 23] and it 
was generally proven in [24] that they are actually equivalent to the radial evolution of the 
classical equation of motion [17]. 

Since the Wilson renormalization group flow theory does not require an ultraviolet com- 
pletion of quantum field theory, the authors of [17] also do not insist on an asymptotically 
AdS region, instead they introduce a finite cutoff r c outside the horizon in a general class of 
p + 2-dimensional black hole geometries. The dispersion relation of the gravitational fluc- 
tuations confined inside the cutoff is shown at long wavelengths to be that of a linearized 
p + 1-dimensional Navier-Stokes (NS) fluid living on the cutoff surface [17, 25, 26]. This 
remarkable relation was investigated in [25] that a given solution of the incompressible 
NS equations maps to a unique solution of the vacuum Einstein equations. An algorithm 
was presented in [26] for systematically reconstructing a solution for the p + 2-dimensional 
vacuum Einstein equations from a p + 1-dimensional fluid, to arbitrary order by extending 
the non-relativistic hydrodynamic expansion proposed in [25]. 

Clearly it is of great interest to develop a holography fluid description dual to an 
asymptotically flat gravitational configuration by introducing a finite cutoff. In this pa- 
per, however, we discuss this issue in asymptotically AdS gravity by introducing a finite 
radial cutoff, because when one takes the cutoff to be infinity, the dual field theory on the 
AdS boundary is well-defined and some results are comparable to those in the literatures. 
Clearly such a study could be a service to further discuss the case in asymptotically flat 
spacetimes. 

Note that one can construct the stress energy tensor of the dual fluid order by order from 
the bulk gravity solution [27]. In this paper we follow the procedures in [14, 17, 25, 26], by 
introducing a finite cutoff surface S c outside black hole horizon, generally discuss the effect 
of finite perturbations of the extrinsic curvature of S c while keeping the intrinsic metric 
of the cutoff surface flat [25]. By applying two finite diffeomorphism transformations, in 
the non-relativistic hydrodynamic expansion limit we obtain black brane solutions, up to 
second order, of the non-relativistic hydrodynamic expansion parameter e, between the 
cutoff surface and the horizon in Einstein gravity with a negative cosmological constant. 
We calculate the stress energy tensor of the fluid on the cutoff surface. The results show 
that the ratio rj/s is still 1/47T, independent of the cutoff, which implies that it does not 
run along the radial coordinate. And it turns out that the stress energy tensor of the 
fluid obeys the incompressible Navier-Stokes equations. We also discuss the case of Gauss- 
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Bonnet gravity with a negative cosmological constant. 

The paper is organized as follows. In Sec. (2) we introduce two finite diffeomorphism 
transformations to a general metric dual to fluid in flat spacetime, while keeping the in- 
duced metric of the cutoff surface invariant. We make the non-relativistic hydrodynamic 
expansion and solve gravitational equations to the second order of the expansion parameter. 
In Sec. (3) we apply this formulism to Einstein gravity with a negative cosmological con- 
stant. We find that the ratio of shear viscosity to entropy density of dual fluid on the finite 
cutoff surface rj/s — 1/47T, independent of the cutoff and that the conservation equation of 
the stress energy tensor of the dual fluid gives an incompressible Navier-Stokes equation. 
In Sec. (4) we consider the case of Gauss-Bonnet gravity with a negative cosmological con- 
stant. The ratio of shear viscosity to entropy density is found to be rj/s = (1 — 8a)/47r, 
and corresponding incompressible Navier-Stokes equations are also obtained there. The 
conclusions and some discussions are included in Sec. (5). 



2 Non-relativistic hydrodynamic expansion 

To study the dynamics of fluid in p + 1-dimensional flat spacetime, we consider a generic 
p + 2-dimensional metric: 

ds 2 +2 = —h(r)dr 2 + 2drdr + a(r)dxidx\ (1) 

where h(r) and a(r) are two functions of radial coordinate r. Introducing a finite cutoff 
surface E c at r = r c (outside black hole horizon if the horizon is present), the induced 
metric on the surface is flat 

r ) a bdx a dx b = —h(r c )dr 2 + a(r c )dxidx\ (2) 
where x a ~ (r, x l ). Introduce proper intrinsic coordinates x a ~ (f, x l ) on S c as 

x° = r= x/M^yr, x l = (3) 
the induced metric is simply given by 

ds 2 +1 = r] ab dx a dx b = —df 2 + 5ijdx l dxK (4) 

In order to keep the intrinsic metric of S c flat, following [26] we take two finite diffeomor- 
phism transformations. The first one is a Lorentz boost with a constant boost parameter 
(3i. In the (f, x l ) coordinates, it is given by 

t^tt-tA^, x i ^£ i - 7 ^r+( 7 -l)^F, (5) 

where 

n a = 7 (l, n 7= 7=i= , A = M J (6) 

\/l - P z 
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In the (r, x % ) coordinates, we have 



u a = . v 2 = = a(r c )5 ljV V , w « = ^/^. (7) 



Thus we can get the boosted metric 

2 °^ r2 ( a ( r ) ?\ 

dsf, 2 = r— - — r -Mr) H — - — -v + 2^drdr - — — -dx dr 

-v 2 /h(r c ) \ a[r c ) J h(r c ^ 

2vi f h{r) a(r) 



1 - v 2 /h(r c ) \h(r c ) a(r c ) 
a(r)S, 



ViVj 1 h(r) a(r) 

" h(r c )(l-v 2 /h(r c )) ~ ^) 



dx l dr 

n(r\ \ 1 



The second is a transformation of r and associated re-scalings of r and 



where we consider the case with /c(r) being a linear function of r as /c(r) = 6r + c, with b 
and c two constants. In this case the general metric (1) becomes 



ds 2 = -h[k(r)} ^rr^ dr 2 + 2bJ ^ r . c \. drdr + a[k(r)} ^ c \. d Xl dx\ (10) 
/#( r c)J y h[k[r c )\ a[k{r c )\ 

When the solution describes a black brane, the cutoff surface r c is required to be outside 
the horizon r c > 7> t . When h(r) = r, a(r) = 1, the metric (1) just describes a flat space- 
time written in the ingoing Rindler coordinates, and the transformations (5) and (9) agree 
with those in [26] if we take k(r) = r — r^. 

After taking the two coordinate transformations one after another, the resulted metric 
still solves the corresponding gravitational field equations. But if we further promote Vi 
and Sk(r) = k(r) — r = {b — 1) r + c to be dependent on the coordinates x a , (that is 
Vi,b,c are no longer constants), the transformed metric is no longer an exact solution of 
gravitational field equations. In order to solve the gravity equations of motion, we take 
the so-called hydrodynamics expansion and non-relativistic limit. Namely we will take the 
scaling 

<9 r ~e°, d^e\ d T ~e 2 ij = l,..p (11) 

together with 

Vi ~ e, Sk(r) ~ e 2 , (12) 

where e will be viewed as an expansion parameter. 

As r is arbitrary between and r c , we demand both (6 — 1) and c scale as e 2 . Then 
up to order e 2 , one has 

h[k(r)} = h{r) + h'(r)5k(r), a[k(r)} = a(r) + a'(r)5k(r), (13) 
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and the transformed metric changes to 
dSp +2 = — h(r)dr 2 + 2drdr + a(r)dxidx l 



+ 



a(r) h(r) 
a{r c ) h(r c ) 

a(r) h(r) 

a{r c ) h(r c ) 



Vidx l dr — 2 



h(r c 



-dx l dr 



v 2 dr 2 



- h(r) 
+ a(r) 



h(r c 

h'(r)6k(r) h'(r c )5k(r c ) 



+ 



h(r c 



-drdr 



h(r) h(r c ) 
a'(r)Sk(r) a'(r c )Sk(r c 
a(r) a(r c ) 



dr 2 + 2 (6 - 1) - b 



h'(r c )5k(r c 
2h{r c ) 



drdr 



dxidx 1 + 0(e 3 ). 



(14) 



The first and second lines of metric (14) are of order e° and e 1 respectively, the other lines 
are all of order e 2 . Note that if one takes 



h(r) = r, a(r) = 1 k(r) = r — 2P, (15) 
and the corresponding non-relativistic scaling 

v i = v e i (T,x i ) = ev i (e 2 T,ex i ), P = P e (r, x l ) = e 2 P(e 2 T J ex'), (16) 

it is easy to see that the metric (14) is the same as the one in [25], up to order e 2 . 

Next we consider the general asymptotically AdS black brane solutions. Following [27], 
we take the AdS radius to be unit so that 

h(r) = r 2 fir), a(r) = r 2 , (17) 

where f(r) is an arbitrary function of r, but it will be given in (22) and (49) for Einstein 
gravity and Gauss-Bonnet gravity respectively. Consider k(r) as a following transformation 

k{r) = r(l - P) Sk(r) = -rP, (18) 

where P is a small parameter. As will be shown shortly, in fact, the parameter P multiplied 
by a factor r ^ l ^ C ) * s ^ ne P ressure density of the dual fluid. 

Substituting (17) and (18) in (14), the first line of the metric with different f(r) solves 
the equations of motion for the corresponding gravity with a negative cosmological constant 
exactly. The remainder terms in (14) could be treated as the perturbations of the metric, as 

P) will turn out to be the small parameters of the dual non-relativistic fluid. To be more 
specific, we consider Vi = Vi(x\ r) and P = P(x l , r) depending on the coordinates x a , but 
independent of r, and take the non-relativistic hydrodynamic limit in (16), as well as (11). 
Then the metric (14) only solves corresponding gravity equations with the cosmological 
constant A = —EiE±ll a t order e 1 . In order to solve the equations to the next order, we need 
to add correction terms to the metric [26, 27]. Let's consider the constraint equation first 
[27]. At order e 2 , we find that there is only one nontrivial constraint condition: diV % = 0. 
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This will turn out to be the incompressibility of the dual fluid. With this constraint 
condition, it turns out that at order e 2 , the source terms only have tensor modes and only 
the following tensor correction terms need to be added to the metric: 

~,2 



—F(r)(diVj + djVi)dx l dx 3 



(19) 



where F(r) is chosen to cancel the source terms at order e 2 and to keep regular at the 
horizon. In order to keep the induced metric 7 a & invariant, we also need to choose the 
gauge such that F(r c ) = [26]. Then our final metric up to e 2 is: 

dsi +2 = — r 2 f{r)dr 2 + 2drdr + r 2 dxidx l 



- 2- 



+ 1 



1 - 



fir) 
f(r c , 
f(r] 



Vidx l dr — 2 
v 2 dr 2 + 



dx l dr 



+ r 2 f{r) 



f(r c 

rf'ir) _ r c f'(r c ) 
f(r) f(r c ) 



r 2 f(r, 

dx l dx 3 



ViVj 



r 2 f(r 



r 2 f(r c 



-drdr 



f(r c 



+ -F{r) {d iVj + djVi) dx l dx j + 0{e 



(20) 



where the terms in last three lines are all of order e 2 . 



3 Fluid dual to Einstein Gravity 

Consider Einstein gravity with a cosmological constant A = —SSE^H in p + 2 dimensions, 
the equations of motion are given by 1 

Efiv = ~ 2 9fJ - u ^ ~ — ~ 9flu = ®' 
The Einstein's field equations admit the asymptotically AdSp^ black brane solution as 

d s 2 p+2 = -r 2 f(r)dr 2 + 2drdr + r 2 dxidx\ f(r) = 1 - (22) 

Another class of dynamical solutions we are interested here can also be found in the region 
between the cutoff surface E c and the black brane horizon up to order e 2 . The approach to 
find the solution is described in the previous section, it turns out that the corresponding 
F(r) term in (20) is 2 



F(r) = / dx 



r c 



x p ) x 2 f(x) 



(23) 



^^Here we use {[i, v, ■ ■ ■ } to stand for the bulk spacetime indices. 

2 We have checked the form F(r) for the case of < p < 8 by mathematic calculation, and we expect 
that this form is also valid for higher p [29] . 
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Here the boundary condition that F'(r) is regular at r = has been imposed. Additionally, 
the integral upper bound has been chosen to keep ^ ab invariant and this matches the result 
in [28] when we take the cutoff surface to infinity. With the constraint condition diV 1 = 0, 
the metric (20) solves the Einstein's field equations (21) at order e 2 , i.e., E^ u = 0(e 3 ). In 
what follows, we will consider the p = 3 case as a calculation example. The fixed boundary 
condition is just the invariant induced metric ^ ab on E c 



"f ab dx a dx b = —rlf(r c )dT 2 + r 2 [dx\ + dx\ + <£c 2 ) . 



(24) 



With the gravity solution (20), one can calculate the stress energy tensor of dual fluid. 
The Brown- York stress energy tensor Z ab evaluated at the cutoff hyper-surface E c is [17] 



8ttG 



(-f ab K - K ab + C% 



ah) 



(25) 



where K ab is the extrinsic curvature tensor of S c , K is its trace, and C is an ambiguous con- 
stant. It is obvious that Z^ b = g^7a& only have order e° terms. After some straightforward 
calculations, we obtain the stress energy tensor of the dual fluid as 3 



T ab = + T t 



■ ab 



z 



(2) 



ab 



0(e 



where 



Zf b ] dx a dx h 



Z^dx a dx b 



Z ( ^dx a dx b 



<fc ( 
n rJ' e 



f c + 2C) dr 2 
Vidx % dr 



(3r 2 J c P + v 2 ) dr 2 + 



fc 

1 + r 2 J c F> 



VjVj 



(diVj + djVi)dx l dx 3 . 



2C + 



r c f' c 



2/c 



(26) 



doc j dx 



r 2 c Pb l3 



(27) 



dx t dx 3 



Since the conservation equations of the Brown- York stress energy tensor are just the Gauss- 
Codazzi formulas of Einstein's field equations [33, 17]. This means that the conservation 
equations of the dual fluid stress energy tensor are the corresponding constraint equations 
of gravity equations. With the conservation equations of the stress energy tensor, we have 
the first nontrivial equation at order e 2 , 



d a Z n 



—diV 1 



da* = 0. 



(2c 



This agrees with the constraint equation in the gravity side and it turns out that the dual 
fluid is incompressible. Taking this to be the case, to the next order, we have 



d a Z ai 



fc 



Vcfcvfc 



d T Vi - 



U (1 + r 2 J c F[) 

fc 



d^i+rlfjl 



Vc 



diP + v j djV t 



0. (29) 



5 Here we use the units IGirG = 1 and f c = f(r c ), f' c = f'(r c ), F' c = F'{r c ) 1 ■ ■ ■ for short. 
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This would be the constraint equation of the gravity metric of order e 3 and more correction 
terms need to be added to solve the equations of motion of gravity at this order. 

In the infinity boundary limit r c — > oo, as the cutoff surface S c is intrinsic flat, from 
the surface counter-term approach [30], we have to take C = —3 for our AdS*> case in 
order to remove the divergence in the stress energy tensor. In AdS/CFT, the background 
metric h a b for the fluid stress energy tensor (Tab) is redefined by stripping off the divergent 
conformal factor from the boundary metric [31] 



b ab 



lim 



lab 



r c ->oo r 



hh ab (T bc ) = lim A /=77 a6 T i 



6c- 



(30) 



Our result can also recover the result for dual fluid on the boundary in the non-relativistic 
limit [27, 28], we will see this later. To see clearly the properties of the dual fluid, it is 
more convenient to rewrite (27) in the (f, x l ) coordinates: 



ff b ) dx a dx h 



f ( ^dx a di b 



f$dx a dx b 



-2 



6V/< 

Tcf' c 



2C) dr 2 

(3idx l df 



6V/C + 2C 



rJ' 



Tcf ' c (3P + /3 2 ) df 2 ^ Tcf ' c 



fc V /< 

(l + r 2 c f c F' c ) {dify + dj^dx'dff 



doc f d/Oc 



2/c 



(31) 



dx l dx 3 



In general, the stress energy tensor of relativistic fluid in 4-dimensional Minkowski back- 
ground 7 a fe = rj a b can be written as [28, 34] 



Tab = pU a U b + ph ab ~ %Wab ~ (0h ab , 



(32) 



where 



<ab 



U a U b + lab, CT a b 



3 



-dh a b, 



e = d a u a . 



(33) 



In the non-relativistic limit and with the incompressible condition, up to order e 2 , we have 



T 



P+(p + P)(3 2 , f T 



'(p + p)A 



(p + PWj + pS tJ - 7/(4/3, + djfc). (34) 



Comparing our dual fluid stress energy tensor (31) with this form, we can get from T^ in 
(31) the energy density and pressure of dual fluid at order e° as 



Po = -(6v^+2C), Po = 6^/T c + 2C 



Tcf c 



= Po + Po 



rj' c 



(35) 



Further we can obtain from T afe in (31) the r c dependent dynamic viscosity 



V(r c ) = (1 + r 2 J c F' c ) = lA 



16nGr 3 ' 



(36) 



As difij ~ e 2 , we see the viscosity r] c is of order e°, it is better to compare this with the 
background of the fluid entropy density, which is just the background black brane horizon 
entropy density. The entropy density s c associated with the cutoff surface is [17] 

*MrJ = ^ =* 1 = 1- (37) 

We see that the ratio of shear viscosity to entropy density is independent of the cutoff r c , 
which means that the ratio does not run with the cutoff. Note that for the static background 
solution (22), the Hawking temperature T# of the horizon and the local temperature T c 
on the cutoff surface respectively are 

rp _ [r 2 f(r)]\ _r h Th _ 1 r h 

H ~ \ r=r h — —i J c - — F== - — F== — • 

We see that the following thermodynamic relation still holds 

r f Arf 

"o = T c s c = —3==, (39) 



for the dual fluid on the cutoff surface. In addition, the dimensionless coordinate invariant 
diffusivity D c defined in [17] is found to be 

D c = T c ^ = * = (40) 

U S c 47T 

a universal constant. 

Next let's read off in the (f, x l ) coordinates the energy density and pressure of the dual 
fluid, to order e 2 , 

Pc = p + 3u P, Pc = Po + co (l + T -^j P- (41) 

The Hamiltonian constraint on S c would play a role analogous to that of the equation of 
state for a conventional fluid [26]. For the fluid dual to Einstein gravity with a negative 
cosmology constant A, if we define T a b = T a b — 2C'j a b = 2(j a bK— K a b), the Hamiltonian con- 
straint would turn out to be T 2 —pT a bT ab + 8pA = 0. This constraint provides the relation 
between the energy density and pressure of the dual fluid. If we evaluate the Hamiltonian 
constraint using the Brown- York stress tensor (31), the first nontrivial equation is encoun- 
tered at order e 2 , which is just the incompressible condition dif5 l = 0. Taking this to be the 
case and substituting (34) into the constraint, we can get (p + 2C)(p + 3p — AC) + 72 = 
up to order e 2 . It can be checked that the energy density and pressure read off from the 
Brown- York tensor of the dual fluid satisfy this equation of state, by use of (p c ,Pc) i n (41). 
We can also calculate the trace of the stress energy tensor (31), to order e 2 , 

T c = f ab T b = 4 [3 {f l J 2 + f^ 1 ' 2 ) + 2C] + -%P, (42) 
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which in general does not vanish. However, if take the cutoff surface r c — > oo, and consider 
the conformal invariance of dual field theory on the AdS boundary, namely, T c — > 0, we 
can also recover the counterterm factor C = —3. In addition, we find that 

lim p c r\ = 3r£ (1 + 4P) , lim p c r 4 c = r\ (1 + AP) . (43) 

Note that in (43) should be absorbed by the dual fluid stress energy tensor (T a j,) defined 
in the flat spacetime h a b in (30), the results reproduce the holographic values. 
Following [28], we define the pressure density as 

PeE ^ = (l4) ? (44) 



Po + Po V 2/, 

and the kinematic viscosity as 



ojq \r c J 4 \r h J AirT c 



c 



1 



(45) 



Then we can rewrite the conservation equation (29) in the (r, x l ) coordinates with Minkowski 
metric j]^ and 3- velocity /3; as 

— d a f ai = B T /3i - v c d 2 fr + kP c + ptdjPi = 0. (46) 

This is nothing, but precisely the incompressible Navier-Stokes equation up to order e 3 of 
the dual fluid in flat spacetime. 



4 Fluid dual to Gauss-Bonnet Gravity 

The Einstein-Hilbert action with the Gauss-Bonnet term £ GB = (R 2 - AR^R^ + R^ var R^ aT ), 
and a negative cosmological constant A = — p ^ 1 - > in p + 2 dimensions can be written as 

S = TTTT^ / dP+2x V=9 (R-2A + aC GB ) , (47) 
lb7rG J m 

where a is the Gauss-Bonnet coefficient with the same dimension as square of AdS radius 
£. As Lqb is nontrivial for p > 3, we will consider the p = 3 case and take the unit 
AdS radius I = 1 in what follows. With the corresponding surface terms [38, 40], we can 
obtain the equations of motion for the Gauss-Bonnet gravity by varying the action (47) 
with respect to metric 

Rfw — -^.Rq^lu — Qg^u + ocH^ = , 

2 ! (48) 

Hfiv = ^{R^uktRJ 1 ^ — 2RnpuaR pa — 2R^ a R c \ + RR^u) — -g^u^GB- 
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The Gauss-Bonnet gravity with a negative cosmo logical constant is solved by [35, 36] with 
spherical symmetry. Here we consider the 5-dimensional black brane solution written in 
the Eddington-Finkelstin coordinates [36, 37] 




-r 2 f(r)dr 2 + Idrdr + r 2 (dx\ + dx\ + dx\ 




(49) 



Note that when r — > oo, f{r) — > j- (l — y/1 — 8aj . One can define the effective AdS 
radius i e in the Gauss- Bonnet gravity as 



1 + VI - 8a 



4a 



1 - v / T Tr 8 



a 



(50) 



We find that in the case of the Gauss-Bonnet gravity, the correction term (19) and the 
solution (20) also work with p = 3. But the form F(r) has to be replaced by 



F(r) 



dx 



5 a 



x 3 (l-2a [2f(x) +xf'(x)}) 



x 2 f(x) 



(51) 



With the constraint equations div 1 = 0, we have checked that (20) solves the equations of 
motion of the Gauss-Bonnet gravity equations (48) up to order e 2 . Here when r c — >■ oo, 
the form F(r) gives the result in [37]. 

From [39, 40], we can take the Brown- York stress energy tensor for the Gauss-Bonnet 
gravity as 



1 



- ab 



8ttG 



K^ ab - K ab - 2a 3 J ab - J j, 



lab 



^PacdbK Cd 



ab 



where 



Pabcd = Rabcd + 2Rb[cld]a ~ ^R a [cld]b + Rla[cld]b 



(52) 



(53) 



is an intrinsic tensor associated with the flat induced metric 7^. So in our case P a bcd will 
not make any contribution to the stress energy tensor of the dual fluid. Namely, we have 



Tab 



8nG 



[Kj ab ~ K ab - 2a (3 J ab - J^ab) + Cjab] 



with 



J ab = =-(2KK ac K c b + K cd K cd K ab - 2K ac K cd K db - K 2 K t 



^cd 



ab) 



(54) 



(55) 



By use of Gauss-Codazzi equations, the conservation of stress energy tensor (54) can also 
be deduced from the equations (48) of motion of the Gauss-Bonnet gravity [39]. Through 
some straightforward calculations, the stress energy tensor of the dual fluid to the Gauss- 
Bonnet gravity is found to be 



- ab 



rp(0) _|_ rp 

1 ab ~r 1 > 



(1) 
ab 



T, 



(2) 



ab 



0(e 



(56) 
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where 



T^dx a dx b 



T^dx a dx b = -2 (1 - Aaf c ) ^ Vl dx l dr 

V fc 



-6V7c " + 8af c V/cJ r 2 / c dr 2 
6 y/% + 2C- 8af c ^fJ c + (1 - 4a f c ) 



Tcf c 



rj' c 



T^dx a dx b = (1 - Aaf c ) (3r 2 J c P + W 2 ) ^dr 



c J_2 



(57) 



+ (l-4a/ c ) 



/c 



2f e 



r 2 c PSij 



dx l dx 3 



[1 - 2a (2f c + r c f' c )} l -±-^Mk(d lVj + d j v i )dx i dx j 



r c Vfc 



The corresponding conservation equations for the stress energy tensor at order e 2 are 



Qa = (l-4a/ c )r c /' 



4 rf 



v7c v7c '"c 

which gives the incompressible condition again, and to the next order 



(5* 



d a T ai 



(l-4af c )f' c 

Tcfc\/~fc 



f 



C -d'\ + r 2 / c 1 



2/c 



diP + v ] d^v l 



0. 

(59) 

where r/ c = [1 — 2a (2/ c + r c / c )] (1 + r 2 / c F c '), which will turn out to be the fluid's dynamic 
viscosity. The equation (59) gives the constraint condition of the solutions of the Gauss- 
Bonnet gravity (48) at order e 3 . 

In the (f, x % ) coordinates, the corresponding stress energy tensor of the dual fluid is 



f^dx a dx b = (-6v^ - 2C + 8af cV ^j0j dr 2 



f^dx a dx b 



f^dx a dx b = (1 - Aaf c ) (3P + (5 2 ) '-y^df 



Q^J C + 2C- 8af c ^/Tc + (1 - 4a/ c ) 
-2 (1 - 4af c ) ^frdx'd? 

V Jc 

fcf'c j~2 



dx 7 dx 



fc 



+ (l-4a/ c ) 



^ L 



Aft 



2/c 



P5 



dx l dx 3 



[1 - 2a (2/ c + r c / c )] (l + r 2 / c F^) + ) 



(60) 
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from which we can obtain the energy density and pressure of the dual fluid at order e° 

Po = -6 V / A-2C + 8a/ cV / A, 
Po = Q^/Tc + 2C - 8af c y/Tc + (l - 4af t 



Tcf c 



(61) 



w = Po + Po = (1 - 4a/ c 



A ri 



while to order e 2 , they are 

Pc = Po + 3w P, Pc = Po + w ( 1 + 
And the trace of the stress tensor (60) is 



2/c 



P. 



(62) 



= 4 [3 (/V 2 + /-V 2 ) - 2a/" 3 / 2 + 2C] + - f° - 



P 



(63) 



If one takes the limit r c — > oo, some surface counter-terms [40] are needed. Since the 
intrinsic Riemann tensor associated with 7 a {, vanishes, in our case, we need only to take 
C = — 2+v/1 ~ 8 ^ ; which is also consistent with a vanishing trace T c = when r c — > oo. Thus 
we have 



lim p c rt = U e r\ (1 + 4P) 



lim p c rl = l e r\ (1 + 4P) 



(64) 



Once again, transforming back to the (r, x l ) coordinates in the boundary metric h a b (30), 
the divergent factor r\ could be absorbed. We have checked the result and shown that 
they match with the stress energy tensor of the dual fluid on the boundary [37]. 
In this case, we find that the dynamic viscosity 



[1 - 2a (2/ e + r c f' c )\ (1 + r 2 J c F' c ) 



'1 -8a). 



The black brane entropy density associated with the cutoff surface is given by 



(65) 



J_r| 
AG r? 



s r Air 



la) 



(66) 



This ratio of shear viscosity to entropy density is the same as that of the dual fluid to 
the Gauss-Bonnet gravity on the boundary [9, 10, 11, 12]. Thus, through the hydrody- 
namic expansion method, we find that the ratio for the dual fluid on the cutoff surface is 
independent of the cutoff, and does not run with the cutoff [14]. In addition, we find the 
thermodynamic relation holds and the diffusivity is changed as 



Wo — T c s c , 



Ik 



An 



(1 -8a) 



(67) 
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where the local temperature on the cutoff surface is 

T _ T H _ 1 [r 2 /(r)]' 1 r h 

Finally we give the incompressible Navior-Stokes equations dual to the Gauss-Bonnet grav- 
ity 

— d a f ai = d T f3i - v c d 2 ^ + ~diP c + ftdjPi = 0. (69) 

OJf) 

where 

fh (1 - 8a) „ /, , rj' c 



"<=t=-Mff> p °={ i+ m p - (7o) 

5 Conclusions 

By use of the non-relativistic hydrodynamic expansion method, we have solved the equa- 
tions of motion for Einstein gravity and Gauss-Bonnet gravity with a negative cosmological 
constant, respectively, and obtained black brane metrics in the region between a finite cut- 
off surface E c and the black brane horizon, up to order e 2 . We have calculated the stress 
energy tensor of the dual fluid on the cutoff surface through the Brown- York tensor on the 
surface [17]. And then we have discussed some properties of the dual fluid. It turns out 
that the dual fluid is an incompressible one in both cases, obeys the Navier-Stokes equa- 
tions, but with different kinematic viscosity. In both cases, the ratio of shear viscosity to 
entropy density is independent of the cutoff, it means that the ratio does not run with the 
cutoff surface. When one moves the cutoff surface to spatial infinity, namely takes the limit 
r c — >■ oo, our results can recover those for dual fluids on the boundary of AdS space [27]. In 
the near horizon limit, our results should also be in contact with the membrane paradigm 
in [41], where the horizon dynamics is described by incompressible Navier-Stokes equations 
in the non-relativistic hydrodynamic limit. 

Our results presented in this paper together with those in [17, 25, 26] show that the 
fluids on a finite cutoff surface always obey the incompressible Navier-Stokes equations, for 
dual gravity solutions up to second order of the non-relativistic hydrodynamic expansion 
parameter. This may reveal some insights for the holographic dual to asymptotically flat 
spacetimes. The study of holography on a finite cutoff may also be helpful to understand 
the microscopic origin of gravity. It is shown in [17] that the entropy flow equation along 
the radial coordinate is equivalent to a radial Einstein equation. We have checked that 
this also holds in the Gauss-Bonnet gravity. 

Finally we would like to stress that the black brane solutions obtained in this paper 
are up to second order of the non-relativistic hydrodynamic expansion parameter e. In 
principle, following [26], we can obtain corresponding black brane solutions to an arbitrary 
order of the expansion parameter. With the resulting solutions, by calculating the corre- 
sponding Brown- York tensor in the gravity side, we can obtain the stress energy tensor of 
dual fluids on the cutoff surface, and then obtain other transport coefficients of the dual 
fluids. 
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